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Abstract 

In this paper we show how the nonlocal effective action for gravity, ob- 
tained after integrating out the matter fields, can be used to compute particle 
production and spectra for different space-time metrics. Applying this tech- 
nique to several examples, we find that the perturbative calculation of the ef- 
fective action up to second order in curvatures yields exactly the same results 
for the total number of particles as the Bogolyubov transformations method, 
in the case of masless scalar fields propagating in a Robertson- Walker space- 
time. Using an adiabatic approximation we also obtain the corresponding 
spectra and compare the results with the traditional WKB approximation. 
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1 Introduction 



In recent years the phenomenon of particle creation from classical sources has ex- 
perienced a growing interest, mainly motivated by its numerous applications in cos- 
mology, but also in other areas of physics. In cosmology, it plays a fundamental role 
in the mechanism of reheating after inflation [|lj which is believed to be responsible 
for the creation of almost all the particles that populate the universe today In the 
reheating models, the oscillations of an homogeneous scalar field (inflaton) around 
the minimum of its potential give rise to an explosive creation of a large amount 
of particles. On the other hand, the same methods are applied to the generation 
of primordial density inhomogeneities in the early universe that later on grew to 
create the present galactic structure |J. In addition, the cosmological expansion 
can give rise to the production of an stochastic background of gravitational waves 
|J. Bounds on the density of these waves are very useful to constraint the different 
cosmological models @]. In all these applications, the method which is used for 
the calculation of the rates and spectra of the particles produced is the traditional 
mode-mixing Bogolyubov technique ||. 

On the other hand the notion of effective action (EA) has proved to be a very 
useful tool for the development of the so called phenomenological lagrangians. Typ- 
ically, effective actions are obtained in theories with heavy and light fields by func- 
tional integration of the heavy modes to find the effective low-energy theory for the 
light modes after some momentum expansion. Usual applications of those techniques 
include low-energy hadron dynamics (the so called Chiral Perturbation Theory), the 
symmetry breaking sector of the standard model, and low-energy quantum gravity 
(see || for a recent review and references therein). Effective actions use to have a 
real and, in general divergent part, that give rise to modifications of the classical 
equations of motion due to quantum effects. Eventually, the corresponding vacuum 
solutions could not exhibit some of the symmetries of the classical theory, thus giving 
rise to the well-known phenomenon of spontaneous symmetry breaking. In addition, 
nonlocal finite terms also appear in the EA which contribute to the imaginary part. 
This imaginary part is physically important since it is connected with the possibility 
of having particle production [[?], |J. By this we mean the production of the quanta 
corresponding to the fields that have been integrated out. 

In this paper, we consider the production of scalar particles from classical gravi- 
tational backgrounds from the effective action point of view. We show how a pertur- 
bative calculation up to second order in the curvatures in the case of masless scalar 
fields, reproduces the well-known general results of particle production in Robertson- 
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Walker space-times and can give rise to the exact amount of particles at least in 
the models we have considered. The paper is organized as follows: in Section 2, we 
review the Euler-Heisenberg lagrangian for QED, but paying special attention to its 
nonlocal part. We show how the perturbative calculation up to second order in the 
coupling constant yields the correct expression for the imaginary part in the masless 
case. In section 3, we introduce the nonlocal gravitational effective action for scalar 
fields and discuss some of the conditions for its application. Section 4 is devoted to 
the actual calculation of the total number of particles produced due to the expan- 
sion in several Robertson- Walker models and the results are compared with those 
obtained by the Bogolyubov technique. In Section 5 we study how to obtain the 
spectrum of the particles and compare the results with the WKB approximation. 
Finally, Section 6 contains the main conclusions of the work. 



2 The nonlocal Euler-Heisenberg lagrangian 

Let us consider the well-know Euler-Heisenberg lagrangian for QED in flat space- 
time || . When the momentum p of photons is much smaller than the electron mass 
M, the one-loop effects, such as vacuum polarization, can be taken into account by 
adding local non-linear terms to the classical electromagnetic lagrangian. Consider 
the QED effective action given by: 

e iW[A] = J [#] [#]e"3 / dixF ^ F " v exp (i J d A xi){% p-M + ie)ip) 

= e-zf^^detiifl-M + ie) (1) 
where as usual p = 7 M (<9 A1 — ieA^). From ([!]) we can write the effective action as: 

W[A]=--J d 4 xF^F^ - z'Trlog((i fl-M + ie)) (2) 
Expanding in a formal way the logarithm we obtain: 

W[A] = — j d 4 xF^ + i ^-Tr[{i? ~ M)~ l Af (3) 
4 J k=i fc 

Using dimensional regularization, it is possible to find the following expression up 
to quadratic terms in the photon field: 



W[A] = J d*x 



1 e 2 2e 2 / 2M 2 

4 ^ 3(4tt)2 ^ (Any { 3 □ 
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- ~fl-2^F(-D;M 2 )W +0{A') (4) 

where A = N e — log(M 2 //i 2 ), N t = 2/e — 7 + log47r is the well known constant 
appearing in dimensional regularization and we have used the expression: 

F{-n- M 2 )F, v (x) = J d*y^e^-^F(p 2 ; M 2 )F, v (y) (5) 



,1 / „2 

,2. j/2'i _ <i / — , / 1 ^ 



with: 

F(p 2 ; M 2 ) = 2 + J q dt log ^ - ^t(l - t) ) (6) 
In the p 2 > 4M 2 case, this function can be written as: 



^ 2 ;M 2 ) = Ji-^-io g V A - (7) 



P 2 C \/T3ip£_| 



In a similar way, the inverse operator 1/D can be defined with the usual boundary 
conditions on the fields as: 

1 -F, u (x) = [ d%^e^-^—K u (y) (8) 



-□ ^ v ' J *(2vr) 4 p 2 + ie 

The expression (|j) for the EA has a regular massless limit. In fact, for small p 
compared with M, the Mandelstam function F(p 2 ; M 2 ) behaves as: 

/ m 2 \ 

F(p 2 ; M 2 ) = - log + 0(M 2 ) (9) 

V — p z — it) 

From (U) we can see that the only contributions in the massless limit are those 
coming, on one hand from the A factor and, on the other hand, from the Mandelstam 
function. Both logarithmic contributions equal, up to sign, so that they cancel each 
other and we obtain: 

W[A) = J d*x (-\F» V F» V - -^-F^T^fJ) + 0(A 4 ) (10) 
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where we have used the following notation: 




T(n)=N f -\oR — (11) 



to be understood as in the previous cases through the corresponding Fourier trans- 
form, with the ie factor as shown in (P). We see that the masless limit of the EA is 
a nonlocal but analytical functional in the gauge curvatures F^ v . 

The EA (U) allows us to derive in an exact fashion the photon two-point one 
loop Green functions. This, in turn, allows us to obtain for example the vacuum 
polarization. In the massive case, the EA can be expanded as a power series in 
p 2 /M 2 , and also in A to obtain the well-known local Euler-Heisenberg lagrangian 



1 e 2 e 2 

£ ,, = f^F AF^F F^UF 

eff 4 ^ 3(4tt) 2 ^ 15(Att) 2 M 2 ^ 

+ 90(4^ ((^) a + j(^) a )+^(^) 8 + ^ 6 ) ^ 

The EA (|J) possesses a non-vanishing imaginary part coming from the Mandel- 
stam function (|7|). This imaginary part provides the pair production rate ||. In the 



massless case (jTOQ we get: 



2 

Im W[A] = 1m J d 4 xC eff = -±- J d 4 x d 4 y—^-e^-^F^(x)F^(y)9( P 2 ) (13) 



(2vr) 



where: 



x < 

0(z) = { 1/2 x = (14) 

1 x > 

The 1/2 value arises as a consequence of the — ie factor in For constant electric 
fields and in absence of magnetic fields, the previous expression gives the probability 
per unit time and unit volume that at least one electron-positron pair is created by 
the electric field: 

p * 21m C eff = ^E 2 (15) 
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Let us compare this result with exact expression for the imaginary part obtained by 
Schwinger ||: 

2 771 2 OO 1 2 
C -£_/ ^ — ^ -L m 717T . . 

P - 2Im ^e// = -^g £B ( 16 ) 

The dependence in the electric fields appears in both a quadratic term and an non- 
analytical contribution exp(— m 2 n7i/eE). This latter term shows the importance of 
the non-perturbative effects in the particle production phenomenon [I]]. However, in 
the masless limit the non-analytical pieces disappear and the result exactly agrees 
with the perturbative calculation in flI5|). Notice that in this case, gauge invariance 
and the dimension of the effective lagrangian constraint the result to be quadratic in 
eE and that is the reason why the second order perturbative calculation gives rise to 
the exact result. Accordingly, in the masless limit, the perturbative calculation can 
provide, in some cases, all the relevant information about the particle production 
processes. In the gravitational case that we will study in the next sections, we will 
show that the same effect takes place. 



3 The effective action for gravity 

Let us consider a real scalar field in a curved space-time with an arbitrary non- 
minimal coupling to the curvature. The corresponding classical action is given by: 

S[</>] = -\j d'x^ (□ + m 2 + £R) <j> (17) 

where: 

□0 = g^V,d u cf> = -Lfy (g^y/gdrf) (18) 

The EA for the gravitational fields that arises after integrating out the real scalar 
matter fields is given by the following expression in Lorentzian signature: 

(0,out|0,in) = Z[g^} = = j [d(j)]e iS[9 ^ 

= J[d(j)}e-^ dx ^ n+m2+ ^ R - ie ^ = (detO)"3 (19) 

where O xy (m 2 ) = {—O y — m 2 — £R(y) + ie)5°(x, y) with S°(x, y) being the covariant 
delta S°(x, y) = g^ 1 ^ 2 (x)5(x,y). Thus we see that, following the analogy with flat 
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space-time we could interpret Z[g^] as the vacuum persistence amplitude. Thus we 
have 

W\g^\ = l -\og det 0(m 2 ) = '-Tt log 0(m 2 ) (20) 

In this expression we have integrated the scalars but the gravitational field is treated 
classically. Accordingly, this EA must be added to the classical action for the grav- 
itational field and it includes the quantum effects due to the matter fields. In 
addition, (|19|) is the generating functional of the Green functions containing scalar 
loops only and external gravitational legs. 

Once one knows the EA at least in some limit, we have all the information 
concerning the semi-classical gravitational evolution in this limit. As we mentioned 
in the introduction, the EA could have a non- vanishing imaginary part, which is 
related to the pair production probability. In fact, the probability P that at least 
one pair particle-antiparticle is created by the gravitational field is given by ||: 

P = 1 - |(0,owt|0,m}|^ = 1 - \ e iW[9 ^ ] \ 2 = l- e - 2ImW ^ (21) 

for small values of Wf^] we have: 

P ~ 21m W[g^\ (22) 

Concerning the applicability of this equation, let us compare the EA method 
with the traditional Bogolyubov technique. The classical equations of motion for 
the scalar field are: 

(□ + m 2 + (R)<j> = (23) 

Unlike flat space-time, there is no natural set of mode solutions to this equation, 
rather we can expand its solutions in different ways, i.e: 

= S(ofcWjfe + 4 w fc) = J2(®kU k + 4«a) ( 24 ) 
k k 

Each of these expansions will give rise to different Fock spaces when interpreting 
the coefficients a*, a[ and creation and annhilation operators. A problem 

arises when we try to identify which of these Fock spaces corresponds to our usual 
notion of particle. In general, this question can only be answered when we have a 
high degree of symmetry (conformal invariance) or if the space-time is flat in the 
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asymptotic in and out regions. However in most of the interesting situations, these 
two conditions are not present. A solution to this problem was suggested in a series 
of works (see |J and references therein) in which the notion of adiabatic vacuum is 
introduced. In the cosmological space-times in which we will be mainly interested, 
in order to define an adiabatic vacuum it is only required that asymptotically in 
the past and in the future the rate of expansion vanishes, i.e. a/a — > with a(t) 
the universe scale factor. Expressing this statement in a covariant way, it would be 
equivalent to require that the curvatures and all their covariant derivatives vanish 
in the far past and future. 

In the effective action approach, Zlg^] can be interpreted as vacuum persis- 
tence amplitude in principle only when the vacuum states \0,in) and \0,out) can 
be defined in regions with a temporal separation [HJ. When this does not occur, it 



is not obvious what is the interpretation of the effective action. However, we will 
show in the following, that the naive calculation of the effective action, in those 
situations in which an adiabatic vaccum can be defined although the space-time is 
not asymptotically Minkowskian, yield the same result for the particle production 
as the standard Bogolyubov technique. As a consequence, in these cases, we could 
try to interpret Z[g^ u ] as adiabatic vacuum persistence amplitude. 

The nonlocal effective action for gravity has been evaluated in different works 
using several techniques. Thus in [JTTJ] it was suggested what would be the form of 



the two-point form factors. In [13] the effective action is derived by means of the so 



called covariant perturbation theory, valid in asymptotically flat manifolds, in JL3 
the same result is obtained by means of the partial resummation of the Schwinger- 
DeWitt series. The result in all these cases up to second order in curvatures can be 
written in the masless case as: 

W[g, v ] = ^/rf 4 x v ^(^^(x)r(D)^ v (x)-^^(x)r(D)^(x) 

+ l(l-t) 2 R(x)T(n)R(x^+0(K 3 ) (25) 

where the form factor T(D) is given in (|TT]) . The local finite pieces as usual depend on 
the different renormalization squemes and they are not relevant for our calculations, 
although in general their coefficients are important to fix the form of the linear terms 
in the trace anomaly. The nonlocal contributions are in any case unambiguous. An 
appropriate representation of the nonlocal form factors is provided by the use of the 
Riemann normal coordinates (the details of this approach will be given elsewhere 



14]). Thus, taking normal coordinates (x^) with origin at yo the action of the form 



S 



factors is understood through the correspondig Fourier transform: 

K(y ) log (^\ K(y ) = J d 4 x^e^n(y Q ) log {^^^j K{x) (26) 

and TZ denotes generically the scalar curvature, the Ricci or Riemann tensors. 

For the sake of simplicity we will study massless scalar particles propagating in 
a cosmological background, whose metric is that of Friedmann-Robertson- Walker 
(FRW): 

ds * = dt 2 _ a{t) 2 ( ^ Kt2 + r 2 d6 2 + r 2 sin 2 (fl)# 2 ) (27) 

where K determines the spatial curvature sign [|15j and a(t) is the universe scale 
factor. 

The EA imaginary part comes from the logarithms in (f25|). Due to the homo- 
geneity and isotropy of space in the present case, the different curvatures appearing 
in that expression only depend on the time coordinate. Thus, we can perform the 
spatial coordinates integration in (|26"D and generically we will obtain: 

Im / d*x-0- 4 e^K(y o ) log {^^) K<?) 

= Im/ dJ^Le^Kivo) log (^^-) 

-TrK{y )K{y ) (28) 

Let us momentarily consider a general metric, not necessarily FRW. It is easy to 
see from the first term in this equation that when the metric is static, i.e only 
depending on spatial coordinates, the argument in the logarithm would only contain 
p 2 — ie. Therefore the imaginary part would be zero and we would recover the well- 
known result of absence of particle production in general (inhomogeneous) static 
backgrounds. 

• FRW metrics 

Returning to the FRW metric we obtain from ( p8| ) the general expression: 

Im W\gJ[ = — [ d 4 iV9 ( —R^ xp Ruu\ P - —R^Rnu 
[ypl 32nJ vy U80 M p 180 " 

2 



1 /l 



+ R 2 \+0(K A ) (29) 
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This result is only valid for homogeneous and isotropic metrics. Comparing 
this result with the divergences, we see that both have the same form. Notice 
that ( p9|) is a linear combination of R 2 u xp^ R\ v and R 2 , but we can choose a 
different basis to write it. In particular, we can take the one made out of R 2 , 
C 2 uXp and E, where C^ u \ p is the Weyl tensor and E = R^xp ~ ^Rfiv + R 2 * s 
the Gauss-Bonnet term. In this basis we have: 

a x R\ v x p + a 2 R\ v + a 3 R 2 = -(ax + ^a 2 ) E + (2a x + C 2 

+ (-ai + ^a 2 + a 3 ) R 2 (30) 

In our case, ai = —a 2 . On the other hand, the FRW metric is locally confor- 
mal to the Minkowski metric and hence its Weyl tensor vanishes. Therefore 
(|29|) only contains the scalar curvature and the Gauss-Bonnet terms, but the 
integral of the latter also vanishes in the class of asymptotically flat metrics. 
Moreover, the asymptotic flatness is not a neccessary condition for the Gauss- 
Bonnet term to vanish and, in fact, examples can be found which are not 
asymptotically flat, but still they have a zero Gauss-Bonnet term contribution 
(see below). To summarize, the imaginary part in these cases reduces to: 

ImWM = ^/^V^Q-e)' l R 2 + 0(n 3 ) (31) 

• Conformal coupling 

In the conformal case (^ = 1/6) it is evident from the above expression that the 
EA imaginary part is zero and accordingly there will be no particle production. 
This is a well-known result and has been proved by studying the positive- 
energy modes of the corresponding Klein-Gordon equation [16| for the scalar 
field. The EA provides in this simple way to prove a general result. 

• Radiation dominated universe 

But conformal invariance is not the only case in which there is no particle 
production in a FRW background. From the above arguments, we have seen 
that the only piece contributing to the EA imaginary part is the R 2 term. If 
this term vanishes, there would not be particle creation. For a FRW metric 
with K = this implies the following condition: 

H = -2H 2 (32) 
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where H = a/a is the Hubble parameter. The solution is simply: a(t) = A(t — 
to) 1 / 2 with A and to arbitrary constants. In fact taking traces in the Einstein 
equations (with the stress tensor corresponding to a perfect fluid), it is obvious 
that R = implies R = 8nG(3p—p) = 0, with p and p the pressure and density 
of the fluid. Accordingly p = 3p, which is nothing but the state equation for a 
fluid of highly relativistic particles. Therefore a radiation dominated universe 
is a stable solution of Einstein equations against pair emission. This result was 



obtained in |16], [17], [u| by means of the Bogolyubov technique, where in order 
to circumvent the problem of the initial singularity, it was assumed that when 
t — > tp with tp the Planck time, the scale factor smoothly tends to a constant. 
This allows us to define an initial vacuum state in the problem. Again the 
out vacuum is chosen as an adiabatic vacuum. Notice that in this case the 
Gauss-Bonnet term can also be neglected. 

Homogeneous anisotropic metrics 

Consider now a general homogeneous but anisotropic metric of the Bianchi 
type I: 

ds 2 = C 2 (r])(dr] 2 - gij { V )dx l dx j ) (33) 

where the 3-metric gij only depends on the time coordinate. Since, as it 
happened with the FRW, the curvatures only depend on the time coordinate, 
it is possible to explicitly perform the spatial coordinate integration in (|26|). 
Therefore we obtain the same combination of curvature tensors as in (|29|) 
for the EA imaginary part. In this case the metric is not conformal to the 
Minkowski one and accordingly it is not possible to drop the Weyl term from 
(|3"0|). The Gauss-Bonnet term continues vanishing under the same assumptions 
about the metric. To summarize, the resulting EA imaginary part can be 
written for this kind of metrics as: 

= ^ / d'x^g [^P^C^ 
+ 0{TZ 3 ) (34) 




This result agrees with that of Zel'dovich and Starobinski [|l9j (see also pOfl ) 
obtained by using standard Bogolyubov techniques. In fact, assuming g^ = 
5ij(l+hi(rj)), neglecting terms of order 0(h 3 ) in fl34]) and imposing that asymp- 
totically the anisotropies vanish, we recover their results. 
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From the above expressions we can extract another consequence. Particle pro- 
duction only takes place when curvature is non-vanishing, i.e, in the presence of 
a genuine gravitational field and not merely by means of a coordinate change as 
it happens for an accelerated observer ||, in the latter case the creation could be 
considered as fictitious. Therefore, for the boundary conditions in the space-time 
geometry that we mentioned before, the EA provides an invariant criterium (inde- 
pendent of the observer) to decide when particle production takes place. 



4 Specific examples with minimal coupling 

In order to illustrate the previous results we will show several examples in which 
the EA allows us to make physical predictions. In some cases it will be possible to 
compare these results with those obtained by means of the traditional Bogolyubov 
transformations. Exact results from the Bogolyubov transformation have been ob- 
tained for a very limited number of models in the literature. 

• Model 1 

We will now consider a complex scalar field and the FRW metric with K = 0. 
It will be useful, in order to compare with other results, to work with the new time 
coordinate defined by: 

T= f a- 3 {t')dt' (35) 

that allows us to write the d'Alembertian operator acting on time dependent func- 
tions as: 

□/(t) = \d T d T f(r) (36) 



First we consider the model proposed in [18|. The scale factor is given by: 



a 4 (r) ~ a\ + e r/s ([a\ - af)(e r/s + l) + b) (e T ' s + 1)~ 2 (37) 

For r — ► oo(— oo), o(r) smoothly tends to a constant a 2 (ai), i.e, it is possible to un- 
ambiguously define initial and final vacuum states. On the other hand, a±, ci2, s and 
b are arbitrary parameters. For a 2 >> ai and using quantum mechanics methods, 
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it is possible to calculate the Bogolyubov coefficients and hence the number density 
of produced particles |18j: 

< N k >= eA J ik _ i (38) 

where the number of created pairs per unit coordinate volume and unit trimomentum 
volume in the k mode is related to the Bogolyubov coefficients by means of: 
< N k >= \f3 k \ 2 § 

The relation between < N k > and the pair production probability per unit 
coordinate volume Pbog is given by this expression [ f2~l|j : 

f d h 

Pbog ^ J [± log(l± < N k >)] (39) 

where the + sign is used for bosons and — for fermions. 

In the present model < N k > does not depend on nor b. Using ( p9|) we find 
for the probability density: 

»°° = / (S l0E(1+ < N - >} = J dk & los ( l + ?=br) (40) 

On the other hand, the EA method provides from (pT|): 



36 32tt 

with: 



p EA ~ 2^-^ / dra 6 (T)R 2 (r) (41) 



a 2 „ a 



R(t) = -12— + 6— (42) 
a s a" 



It is possible to perform the integrals in ((Sj) and (|41| ) in an explicit way so that 
we can compare both results at the analytical level. They yield exactly the same 
result: 

PEA = PBOG = 5760^p (43) 
where we have taken the limits 6 = and 02 - > 00. 

• Model 2 
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The second model we will study is that proposed in |22 
given by: 

a\r) = A 2 r 2 + B 2 



The scale factor is now 
(44) 



where A and B are arbitrary constants. In this case, space-time is not asymptotically 
flat and therefore the Bogolyubov calculation is based on the definition of adiabatic 
vacua. However, the Gauss-Bonnet contribution vanishes and thus we can use again 
The number of created pairs in the k mode is given by: 

< N k >= e — (45) 
Once again both methods yield the same results for the probability densities: 

7v4 3 

PBGC = PEA = (46) 

Since in ([IT]) we have neglected higher order terms in curvatures, we can conclude 
that in these two cases they do not contribute to the EA imaginary part. As we 
found in the QED case, here again the second order perturbative calculation is exact. 
To check this fact, we should calculate the complete expression for the EA as we 
did in Section 2, however the very same arguments used in that section suggest 
that in the absence of a mass term, since in both cases there is just one dimensional 
parameter, it is not possible to build any other term with the appropriate dimension. 



5 Spectrum and WKB approximation 

The traditional Bogolyubov method for particle production gives information, not 
only on the total number of created particles, but also on their energy distribution. 
However, only in very specific cases, closed analytical expressions can be written. 
As we have seen, the EA method provides a closed expression for the total number 
of particle that is obtained from the curvatures and, therefore, can be evaluated for 
arbitrary scale factors in a very easy way. In this respect the EA method is obviously 
more advantageous than the Bogolyubov method. However, it is not obvious how 
to derive the spectra in this formalism. Let us try to clarify this issue with a simple 
example and compare our result with the one obtained from the traditional method. 

Consider the Klein-Gordon equation for a minimally coupled massless complex 
scalar field 

□0 = (47) 
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Introducing the FRW metric fl2"T|) with time coordinate r and K = 0, we look for 
solutions by means of variable separation <p(r,x) = Xk{T~)e lkx . Hence the temporal 
equation can be written as: 

d 2 Xk 



dr 7 



+ a\r)k 2 X k = (4£ 



where k 2 = kk. In the simple example we are going to consider, the scale factor is 
made of two step functions: 

G 4 (r) = 1 + v 2 (6(t + T) -6{t-T)) (49) 

with v and T being arbitrary parameters. The Bogolyubov coefficients provide the 
following value for the number of created pairs per unit coordinate volume and unit 
trimomentum volume in the k mode: 

< N k >= \/3 k \ 2 = 4(i ^ 2) sin 2 (2TkVT+^) (50) 
Expanding the RHS of (]39D using ( |50f) up to 0(v 4 ) we find: 



p BOG ~ dkk 2 sm 2 (2kT) (51) 

£>7T Z JO 

The integrand gives the probability density per unit trimomentum volume. On the 
other hand, the EA method gives the following result from (|3TD : 

1 1 r 

Pea ^ 2Imw ~ 4—— / dra 6 (r)R 2 

= "wAI^-v-^^ (52) 

We have introduced a global 2 factor in the EA because now the field is complex. 
The spectrum can be obtained by introducing a complete set of plane waves: 



p EA ~ 4—— / drdr' f°° ^a 3 (T)R(T)a 3 (T')R(T')e' 
y 32tt 72 J J-oo2n K ' K J K ' K ' 



■ipfr—T 1 ) 



V 



1 



647T 2 Jo 



oo 



cipp sin (pT) (53) 



Comparing pbog with we find that both integrands agree by identifying 
p = 2k. This is sensible and represents the energy conservation in the pair creation, 
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since k is the single particle energy and p is the energy of the gravitational field 
oscillations producing particles. 

From the complementary point of view, given the number density of created 
particles < Nk >, it is also possible to reconstruct the scale factor evolution by 
inverting the previous steps: 

2 

+ 0{v 6 ) 



a 6 (r)R 2 (T) 



36 



dp 



2tt 



Wlog(l+ < N p/2 >)e 



-ipr 



dp sin(pT) 



oo ZTX 



-ipr 



+ 0(v 6 ) 



Q7; 4 

^L(S'(T-T)-8'(T + T)f + 0(v 6 ) 



(54) 

This result agrees with the calculation from (f49|). 

Let us try to generalize the above results for arbitrary scale factor evolution. In 
the above example, it can be shown that the difference in the results using plane 
waves or a complete set of solutions of the Klein-Gordon equation is 0(v G ). Therefore 
the former is a good approximation. Now we have to take into account the presence 
of the curvature. With that purpose we take a variable-frequency plane waves such 
that in the vanishing curvature limit they tend to the usual plane waves. 

Let us consider the temporal part of the Klein-Gordon equation fl4*8|). This is a 



harmonic-oscillator equation but with a time-dependent frequency a>fc(r) = ka (t). 
Changing to the new time coordinate drj = a 2 (r)dr, the equation can be written as: 



d 2 Xk 2 ddxk 



k 2 Xk = 



(55) 



drf a drj 

In the limit in which the expansion rate a /a is much smaller than the frequency of 
the oscillations k, the equation reduces to the flat space-time form. Therefore let us 
consider that limit and let us introduce a complete set of plane waves corresponding 
to the new time coordinate 77: 

„ 1 1 



Pea 



1 1 
) 

'36 32tt 

1 1 



dra b (r)R 2 



36 32vr 



drjdrj' a 2 (r])a 2 (7/) R(r)) R(rf) 



36 32tt 

Changing again to the old coordinate r we have: 



d^^Ri^f 

dp -ip(r,-n') 
2tt 



Pea 



1 



57T JO 



dk 



dr 



-2ik J^a 2 (r')dT' 



(56) 



(57) 
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where we have used (f42|). According to the above discussion, the introduction of the 
plane waves only makes sense in the adiabatic limit and therefore this expression 
is valid only in the limit a/ a << k. As in the step function example, we have 
used p = 2k. The pairs density < Nf~ > can be calculated in a very easy way by 
identifying ( |57|) with fl39|). Obviously from the equality of the integrals we cannot 
obtain the equality of the integrands, however the covariance and dimensionality 
of the integrands allows us to constrain them. In fact, from the equality of (|57|) 
and (|39| ) we know that the integrands differ at most in a function f(p) such that 
J dpfip) — 0- Now let us assume that such function can be written as (up to OilZ 2 )): 

f(p) = f dr ) dr ] 'e ip(T '- T '' ) a 2 {r])a 2 (ri')F(r])F(r]') (58) 



where F(rj) and F(rj') are some appropriate functions. The condition / dpf(p) = 
implies: 

' dr]a 4 (ri)F(r])F(ri) = (59) 

Since the integrand has to be a dimension 4 operator and a scalar function then the 
only possibility satisfying that condition is: 

FF = aE (60) 

where E is the Gauss-Bonnet term defined before, which is a total derivative, and a is 
some arbitrary constant. However we know that in a radiation dominated universe, 
where R = 0, the spectrum is identically zero < N k >= (see |16|, [H], |18|). This 
implies that the contribution from the Gauss-Bonnet term should also vanish. This 



fact allows us to fix the constant a = 0. As a consequence the result in (|57|) gives 
the correct spectrum up to 0{1Z 2 ) at least in the adiabatic limit we are considering. 
In fact, as shown in Table 1 the results are in good agreement with the Bogolyubov 
method specially for large values of k. 

Since we have used an adiabatic approximation in the last step, we can try to 
find which are the differences with respect to the usual WKB approximation in |23| . 
In this method, the solutions of the equation (^) are taken to be: 

Xk (r) = a(r)e~ 1 ^ + P(r)e 1 ^ (61) 

with 

x/j(r) = k f T dr'a 2 {T') (62) 
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with boundary conditons a(— oo) = 1, /3(oo) = 0. Putting this ansatz back into the 
equation of motion, we get: 



ae 



(63) 



where the condition de l ^ + /3e 1 ^ = is used (see |23| for details). The solution to 
first adiabatic order is given by: 



f3{r) = - [°° dr'-e-WW 
Jt a 



The probability density is given by: 



Pwkb 

From (|4]) and (|65[) we see that 

1 



1 /■ d 3 k 



2tt y 4tt 2 



-oo 



(64) 



(65) 



■57T° 
1 '' x 



d 3 k\(3(-oc)\' 



in 3 Jo 



rdk^k 2 [dT~ e -™j; a2 w 

Jo J a 



2vr 
1 



-Jdk [drhe-^f^ 2 ^' 
7r J J a 



2tt 2 



dk 



dr d 1 d c~ 2ik Jo a2{T ' )dT ' 



a —2ia 2 dr 
d ( a 



8tt 2 J J dr W 



2ik F a 2 (r')rfr' 
« 



(66) 



In the last step we have used integration by parts assuming that a /a 3 vanishes for 
t — ► ±oo (which is the same condition as for the vanishing of the Gauss-Bonnet term 
contribution in (|29|)). This condition is satisfied in the models we have considered 
in the paper. Therefore we get: 



87T Jo J \a a a^J 
which is valid again only in the adiabatic limit. 



(67) 
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k(a.u) 


< N k > BOG 


< N k > EA 


< Nk >WKB 


0.3 


5.31 10~ 4 


4.15 10~ 4 


4.96 10" 4 


0.4 


4.31 lO" 5 


3.33 10- 5 


4.13 10~ 5 


0.7 


2.28 lO" 8 


1.75 10~ 8 


2.28 lO" 8 


1.0 


1.22 10~ n 


0.93 HT 11 


1.24 10" 11 


1.2 


7.97 10~ 14 


6.48 10" 14 


8.64 10" 14 


1.3 


6.44 10~ 15 


5.91 HT 15 


7.76 10~ 15 


1.4 


4.44 10~ 16 


4.90 HT 16 


6.19 10~ 16 



Table 1: Number densities corresponding to the model ( p7| ) with s = 2(a.u)~ 1 , 
a\ = 1, a% = 500 and 6 = 0. BOG denotes the Bogolyubov method and EA the 
effective action. 

In Table 1 some values of the number densities are shown for the different meth- 
ods. The results have been obtained from (0) and ( |67| ) for the model (|37| ) by 
numeric integration. Due to the strongly oscillating integrals, the results can only 
be given for small momenta. Both methods give similar results to those obtained 
with the Bogolyubov coefficients. Notice that, despite the fact that the WKB and 
EA expression in ( |5TD and ( |5"7| ) are not identical, there is no contradiction in this. 
These results come from different approximations. The WKB method comes from 
a derivative expansion and is not covariant, whereas the EA is an expansion in cur- 
vatures and covariance is imposed from the beginning. From Table 1, we see that 
the contribution from the (a/ a 2 ) 2 terms in ([57]) and ( |B~7|) is always smaller than the 
contribution from a/a 3 . 

6 Conclusions 

In this work we have shown how to use the nonlocal form of the gravitational 
EA (up to 0(7Z 2 )) for the computation of massless scalar particle production. 
For FRW backgrounds in which the expansion rate asymptotically vanishes, it is 
shown that the particle production probabilities only depend on the scalar cur- 
vature. As a consequence and as expected there is no particle creation in a ra- 
diation dominated universe. This is also the case for conformally coupled scalar 
fields. For anisotropic homogeneous metrics we reobtain the well-known expres- 
sion of Zel'dovich and Starobinski. We compare our results with those obtained by 
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means of the well-known Bogolyubov transformations. In the examples considered, 
the agreement between both methods is complete for the probability densities. This 
fact is quite remarkable since we have used a perturbative expression for the EA, 
whereas the Bogolyubov results are exact. This indicates that in some cases a per- 
turbative calculation may contain all the relevant information about the particle 
production processes. In principle, the EA is defined for asymptotically flat mani- 
folds, however it is interesting to notice that the naive extension to those manifolds 
in which adiabatic vacua can be defined, properly reproduces the correct results. Fi- 
nally we have also compared the different spectra derived with the EA, Bogolyubov 
and WKB techniques. 

In principle the EA method can be extended to more general metrics (not nec- 
essarily homogeneous) in a straightforward way. This fact could make it valuable in 
those areas in which the Bogolyubov technique has been traditionally used. In addi- 
tion this method can also be applied to the production of higher spin particles such 
as Dirac and Weyl fermions, gravitons, gravitinos, etc. Finally, in a recent work [p4 



the relevance of the nonlocal EA for particle creation has also been stressed from a 
different point of view based on the energy-momentum tensor expectation values. 
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